Macroscopic Quantum Tunneling in "Small" Josephson Junctions in magnetic field 
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We study the phenomenon of macroscopic quantum tunneling (MQT) in small Josephson junctions 
(JJ) with an externally applied magnetic field. The latter results in the appearance of the Fraunhofer 
type modulation of the current density along the barrier. The problem of MQT for a point-like JJ 
is reduced to the motion of the quantum particle in the washboard potential. In the case of a finite 
size JJ under consideration, this problem corresponds to a MQT in potential which itself, besides 
the phase, depends on space variables. Finally, the general expression for the crossover temperature 
To between thermally activated and macroscopic quantum tunneling regimes and the escaping time 
t osc have been calculated. 

PACS numbers: 74.50.+r 



The paramount importance of macroscopic quantum tunneling lies in the role of such a process in a variety of 
physical phenomena, which span from the cosmological implications of the false "vacuum" to recent studies on BEC 
systems. Great attention deserves in this context the occurrence of MQT in Josephson junctions^ observed both in low 
T£ and high T<& structures. Usually the phenomenon of MQT is considered in a "small" JJ, i.e. in the assumption 
that the characteristic size L of the junction is much smaller than all other related parameters of the problem (some 
exceptions can be found in experimental* and theoretical papers^ ). Such parameters are the Josephson penetration 

depth Aj = (ftc 2 /87rej 1 ^ off ) 1 ^ 2 and the characteristic length L H = £%/£ e s, with £ H — (he/ eH^) 1 ^ 2 as the standard 
quantum mechanical magnetic length. The external magnetic field H cxt is supposed to be directed along the plane of 
junction. The critical current density j\ of the JJ will be defined explicitly below. The effective length £ e g- depends 
on the relation between the thickness du\ (i = L, R) of the superconductive electrodes and the magnetic penetration 
depth A(,;) of the bulk superconductor materials, which the electrodes are made of. In the limiting cases one can 

f _ f A(i) + A( 2 ) + d xi A(j) < 

cff \ + d(2) + 4x, < A(i) ■ 1 ' 

The finite value of the junction length L becomes important even for L < Aj. Such a situation takes place when the 
magnetic field is applied to the junction what results in appearance of the Fraunhofer diffraction pattern I c vs H ext . 
In our analysis we will suppose that L <C Aj but will not impose any restrictions on the relation between L and Lh- 

The effect of magnetic field on the properties of JJ leads to three different manifestations. First of all magnetic 
applied filed suppresses the value of the order parameter A and changes the frequency dependence of the single- 
particle density of states p (ui) (formally, the presence of magnetic field changes the frequency dependence of the Green 
functions)^ . Then, its presence leads to the spacial modulation of the current flowing through the junction^i^iii. 
Finally, the magnetic field affects on the probability of the electron tunneling through the barrier. Since the thickness 
of insulating barriers in the case of JJ is very small (about 10A) this effect we will neglect in our discussion. 

Let us start our analysis from the suppression of the order parameter A and changes of the frequency dependence 
of the Green functions by the magnetic field. The Green functions integrated over the energy variable ( a, (3) for both 
electrodes can be found from the equation^: 

a 2 

dv 2 




aA-(3u + -{a[—- 2ieA ) (3 - fl—a \ = a@T, (2) 



with the normalization condition a 2 (r,u>) + (3 2 (t,ui) — 1. We choose the vector potential gauge in each electrode in 
the form A = (Hz, 0, 0), where the direction of z axis is supposed to be perpendicular to the junction surface and the 
magnetic field H is directed along the y-axis. The coordinate z in each electrode changes in the symmetric interval 
z G [— d/2, d/2], while the coordinate along the barrier x S [-L/2, L/2]. In view of the smallness of the JJ length one 
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can see that the spacial dependence of A and Green functions can be neglected. Averaging the Eq.Q over the space 
variable one can determine the frequency dependence of the averaged Green functions (a (u>) , (3 (w)) 

5(i)A w - /3 (i) cj = (r (l) + D {l) e 2 H 2 xt d 2 {l) /6) , (3) 

where r^j = t~^\ is the inverse electron spin-flip scattering time. 

The critical current density J\ (H) can be determined from the Ambegaokar-Baratoff type formulapii: 

Ji(H) = -^T^ {1) (cj)0 (2) (-u,), u = u, n = 2nT(n+~\, (4) 

where Rn is the junction resistance in normal state and S is the junction area. Thus the effect of magnetic field results 
just in the renormalization of the critical current density. The latter becomes essential when D^e 2 H 2 xt d 2 ^ ~ A. 

Let us discuss of the second external magnetic field manifestation on the properties of a J J, related to the space 
modulation of current density in it. In the case of a "point-like" J J the value Ji determines the height of the potential 
barriers and hence the escaping time for the MQT. In the case of a small J J of finite size the current density is space 
modulated by the external magnetic field and the escaping time calculation becomes more complicated. Indeed, while 
for the point contact it is reduced to the motion of the quantum particle in the washboard potential in the case under 
consideration the problem corresponds to the tunneling of a lengthy object in the potential which itself, besides the 
phase, depends on space variables. 

To calculate the escaping time r csc we need to write down the functional of effective action A [ip] , what can be done 
in this case basing on the results o f 12 i 13 : 




x a L (r- n) a R (r - n) + cos (ip (r, r) + <p (r, n)) f3 L (r - n) f3 R (r - n)] 

«KI* f 1/2T , sin2 ( y(r ' r) ~/ (r ' Tl) ) h 2 c 2 S (dip(r,r) e4 ffrTJ V 

-f- — S / . i 



2sin 2 ^ (r ' T) ~^ (r ' Tl) 



dT i — . 2 , rr, ^ + z — 2^ — ?r rl H «' xn ] f ■ ( 5 ) 

1/2T sm 2 [itT(t - n)] 8we 2 £ eS \ dr he ' ' 



Here C is the junction capacity, I is the full current flowing through the junction, i? s h is the shunt resistance. 
The integrals are carried out over the imaginary time and the junction area. The first term corresponds to the 
kinetic energy of the junction. The second and fourth terms describe the contribution of the potential energy to the 
action. Let us stress that the third term, corresponding to the capacity renormalization appears only beyond the 
quasi-classical approximation. The action accounting for all these terms was already found i n 12 i 13 in the case of 
microscopic study of the Josephson current decay in point-like junction. The last term in Eq.([5]) accounts for the 
finite size of the junction and the magnetic field contribution to the effective action. The motion in the classically 
permitted region is described by the same functional Eq. ( [5]) in the assumption r = t\ what eliminates the third 
term. The further variation procedure of Eq. ( [5]) over ip results in the getting of the quasi-classical equations of 
motion&l&ii. 

We are looking for the escape time r csc with the accuracy up to the pre-exponential factor. The exponential factor in 
it is determined by the extremal value of the functional Eq. ( [5]) found at the class of solutions periodical in imaginary 
time r with the period l/T. The function ip (x, r) at which the extremum of functional Eq.([5]) is realized (the saddle 
point of functional ) we will search in the form 

x _ 1 „ 

¥> (x, t) = <p - V p (x, t) = -$ (x) + tp (x, t) , (6) 

Lh 2 

where p>o — const, is the phase difference at the junction in the absence of magnetic field, while the function ip (x, r) 
satisfies the boundary condition^ 



dip 
dx 



x=-L/2 



dip 
dx 



= 0. (7) 

x=L/2 



The experimentally realistic case corresponds to the current close to its critical value I c . This is why we restrict 
our consideration namely by the condition I — > I c . Let us note that in this case the third term in Eq.([5]) can be 
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excluded from consideration by the simple renormalization of the junction capacity C 

nh f°° duj ( da.h\ ( dotR 



C*, where 



C* =C 



2i?N 7^oo 27r \ du 



duj 



(8) 



The substitution of Eq.( [6J to the Eq.( [5]) with the further expansion of the cosine over tp (x, t) in the fourth term 
leads to the cancellation of the linear in <Jp second term. Since the function tp has to be periodical in imaginary time 
t with the period 1/T we present it in the form of the Fourier series. This makes the effective action ([5]) to become 
the functional of the Fourier coefficients b n (let us note that b n (x) = (x)): 

hi f 00 

Alb n ] = -±jd 2 r {[cos^(x)} + C(n 2 + Q sh \n\)]bl 



-f sin [$(*)]( £ b k b k+n )b n + X 2 (^- 



\k=-z 



(9) 



where we have introduced the notations C = 2ir 2 T 2 C* / (eHJxS) and the Q-factor Q^ 1 = 2ttT R sh C* / h. The further 
variation of the functional © over b n results in the system of equations: 



d 2 h 



A 2 j^ + [cos [$ (x)} +((n 2 + Q sh \n\)] b n - sin [* (x)} 



(2 - S 7h0 )b b n + 2 ^2 b k br, 



k=l 



E 

fc=i 



bkb„ -I, 



(10) 



Let us recall that we will look for the solution of this system with the boundary conditions ([7]). The analogous system 
of equations for the point junction was solved for the cases of low and high shunt resistance^ and also close to the 
temperature To corresponding to the crossover between the thermally activating (Arrenius law) and quantum regimes 
in tunneling through the barrier. 

In order to calculate the escape time we need to know the value of the effective action ([5j at the saddle point, i.e. 
at the optimal trajectories b n (x) , being solutions of the system ( fTU|) . The related expression can be considerably 
simplified, subtracting from the functional Eq. ( [5]) the Eq. ( ITUl) multiplied by corresponding b n (x) : 



^M = H / A-sin[<JM 



n=l k=l 



n=l k=l 



(11) 



We consider the junction of a small enough size: L C Aj. Thus the system Eq. (|10p acquires the large parameter 
X 2 / L 2 3> 1. Consequently the functions b n (x) weakly vary with the change of its argument and can be represented 
in the form: 



b„ = (b n )+b n , \b n \<^\(b n )\, (12) 

where (...) indicates the average over the coordinate x along the junction. The corresponding equation for b n takes 
the form 

A ^ = (5 " )sin(2 ^ o)sin (li)- (13) 
Its solution satisfying the boundary conditions ( [7]) is 



b n = ttt {b n ) sin (2^ ) 

What concerns the equations for the average values {b n ) they can be obtained directly from the Eq.( fTU|) by the 
method of subsequent approximations. 

Now we can calculate the escape time r csc near the crossover temperature To, taking into account the effects of the 
current density coordinate dependence, due to the applied magnetic field. All coefficients (b n ) with n > 1 are small 
and quickly decrease as n increases. It is why at temperatures close to To in Eq. ( lll[) one can restrict consideration 
by the contribution of 6q, b\ and 62 (~ b 2 ) . 
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The crossover temperature To is characterized by the fact that below it the quantum tunneling through barrier 
becomes more probable than the thermally activated escape process. Formally this temperature is the bifurcation 
point of the system ( fTTJj) when the first non-zero solutions for n > 1 appear. Above it exists the only, trivial, solution 
is bo — const, what results in the action ( [5|) A [ip\ = AU /T corresponding to the Arrenius law for the escape time. 
Below To the non-zero solutions b n appear, what leads to the appearance of the imaginary time dependence of (p , the 
action ( [5]) changes its form and, consequently, the temperature dependence of the escape time r esc ~ exp [—A (<p , L)] 
changes qualitatively. After some cumbersome but straightforward calculations one can find the quadratic equation 
determining the crossover temperature To : 



27r 2 C*T 2 
he 



where I c = Ji(H)S is the value of junction's critical current, while 

2 



T (<Po, L) = cot ( 2< £o) + sin ( 2 <A)) ( 7^- 
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Lh~ 
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(15) 



(16) 



In the case of a point junction the right-hand-side of Eq. ( [151) reduces to \Jl 2 — I 2 . 

In vicinity of To the essential contribution to Eq. ( [TTj) give only the terms with bo, 6i and b^. Their average values 
can be found from Eq.f W2\ -{ fT4^) . We present here only 



(bif = 



T (<A),£)C sin 


(2^ )(^)sin( I 


L ~) C 2 (Qsh + 
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2 sin 2 (2<^o) 




2 1 
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what finally permits to find 



~ r . hJiS . , . ( Lh\ . ( L 



T J {<po,L) — sin (2ip ) I — I sin I — 



(17) 



1 



Q S h + 1 Q S h + 3 



One can see that the effect of magnetic field in a finite size junction, even in the case L -C Aj, does not reduce 
simply to the renormalization of the value of I c . Indeed, both characteristic parameters of the problem, To and the 
escape time r esc , turn out to be the functions of two parameters: I/I c and L/ Lh- Let us stress that when the current 
/ is very close to its critical value I c the approximation (|12p for the functions b n (x) becomes incorrect, and this region 
has to be studied separately. These results are of great interest for the whole scenario of quantum coherence, as for 
instance in the context of MQT in Josephson structures based on d-wave superconductors^ 
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